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Abstract 

In the paper we check whether the contribution of (— )(1 + z) 6 type in the Friedmann equation can be 
tested. We consider some astronomical tests to constrain the density parameters in such models. We describe 
different interpretations of such an additional term: geometric effects of Loop Quantum Cosmology, effects 
of braneworld cosmological models, non-standard cosmological models in metric-affine gravity, and models 
with spinning fluid. Kinematical (or geometrical) tests based on null geodesies are insufficient to separate 
individual matter components when they behave like perfect fluid and scale in the same way. Still, it is 
possible to measure their overall effect. We use recent measurements of the coordinate distances from the 
Fanaroff-Riley type lib (FRIIb) radio galaxy (RG) data, supernovae type la (SNIa) data, baryon oscillation 
peak and cosmic microwave background radiation (CMBR) observations to obtain stronger bounds for the 
contribution of the type considered. We demonstrate that, while p 2 corrections are very small, they can be 
tested by astronomical observations - at least in principle. Bayesian criteria of model selection (the Bayesian 
factor, AIC, and BIC) are used to check if additional parameters are detectable in the present epoch. As 
it turns out, the ACDM model is favoured over the bouncing model driven by loop quantum effects. Or, 
in other words, the bounds obtained from cosmography are very weak, and from the point of view of the 
present data this model is indistinguishable from the ACDM one. 
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I. INTRODUCTION 



Modern astronomy has led us to the cosmological concordance model (CCM). The interpretation 
of observational data suggest that we are living in almost spatially flat, low density, accelerating 
universe filled with dust-like matter . In other words, at first sight there is nothing more 

to the universe than the simple ACDM model. In this model, apart from the observed baryonic 
matter we have two additional components. One is in the form of non-relativistic and cold dust 
matter contributing one third of the total energy density, while the second one, having negative 
pressure violating the strong energy condition, is contributing about two thirds of the total energy 
density. This second term is interpreted as non-zero cosmological constant A. It could be called 
the minimal model in that it is the simplest to fit the data so well. 

Despite that, it is obvious such model is far from reality, and much attention is paid to possible 



extensions of physics in that area 



Especially the early universe, where general relativity 
needs closer association with quantum theories, is open to modelling. In order to preserve the 
agreement with late epoch observations, the modifications need to introduce effects negligible at 
present, and dominating for early times, that is to say, decreasing with the scale factor a, or 
increasing with density p - if one assumes an expanding universe. The Loop Quantum Gravity 
offers the possibility of the description of the discreteness of spacetime at the Planck scale. The 
classical ACDM should be emergent from the more fundamental quantum model. Singularities in 
classical general relativity should be replaced by the quantum description of the very beginning 
of the gravity. For example the Loop Quantum Gravity cosmological model predicts the bounce 
instead of the initial singularity. For full description of the universe evolution both the quantum 
and classical regimes should be taken into account. 

Here, we investigate the observational constraints on a contribution of p 2 type to the Hubble 
parameter H 2 (z), which could also follow from an equivalent term of type (— )(1 + z) 6 . Because we 
have no a priori information which physical theory is valid in the Universe at the Planck epoch all 
possibilities should be taken into account. The viable interpretations of the presence of this term 
in H(a) come from the loop quantum, braneworld and non-Riemannian cosmologies. 

A. Universe in Loop Quantum cosmology 

One of the important, unsolved problem of the modern cosmology is the problem of initial 
conditions for the Universe. Because, at least from our point of view, the Universe is given in one 
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copy the initial conditions cannot be taken from "outside". 

Many cosmologists argue that problem of initial conditions for the Universe can be shifted 
to the Planckian epoch in which quantum effects were crucial. In this context, the idea that 
present expansion of the Universe was preceded by a contracting phase seems to be very attractive. 
Geometric effects in Loop Quantum Cosmology (LQC) predict the presence of p 2 modification to 
the Friedmann equation H 2 (a) of negative sign [g, Q] which modifies the early evolution leading to 
a bounce in the generic case. 

The basic formula H(z) which is used in cosmography assumes very special form in the case 
of loop quantum cosmology if effects of curvature are included. In this case effects of curvature 
cannot be simply included in a model by adding a noninteracting curvature fluid with energy 
density p^ = — ^,Pk = — ^ because loops quantum effects manifest themselves in the form of a 
modification to the classical Friedmann equation. For example if we consider a closed model the 
explicit form is given by [?], 8, 0, 

* = + ^) (i-JL __*_ + _,?_) (1) 

\3 3 a z y V Pcrit Pcrit a * p CT it J 
where p CT it is related with Planck ener gy d ensity p CT n = 16 nf a pp^ — -^ 2L - 



Ppi = is Planck density 
and 7 ~ 0.2375 is the Immirzi parameter. [Ill] . With p CT { t — 0.82p p i ll2| and p p i = 5.155 x lO 96 ^ it 
leads to p cr ;t — 4.19 x 10 93 ^j. Please note that only the first term in parentheses is corresponding 
right hand side of the Friedmann equation. Therefore, loop high-energy modifications of general 
relativity appear in the second term in parentheses. The classical general relativity limit can be 
obtained if p cr [ t goes to infinity and the second term approaches unity. 
The effective Friedmann equation assumes the form 

H 2 = - ( 1 — ) + other terms (2) 

3 V Pcrit/ 1 ' 

The applications of LQG methods to cosmology , offer the possibility of investigating the 

consequences of discreteness of spacetime on the quantum level. One can distinguish a quantum 
bounce from the classical one, for which p <C /o C rit and the quantum corrections are negligible. 

In the next point we pointed out similar p 2 modification of the Friedmann equation as in LQC 
(for the comparison of brane world scenarios and LQG cosmology see 11, 16, 17. Iisl fl^|). 

Recently Bojowald pointed out the possibility of testing loops quantum effects using more 
precise data in the near future [2^. He pointed out that while quantum corrections from gravity 
are commonly expected to be small, a quantum structure of spacetime can produce potentially 
observable effects. Martin Bojowald argued that cosmological models of loops can be in principle 
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testable but there is no concrete proposal in this context. We indicate some possibilities of such 
testing. However, we must remember that the expansion history Htz) measures only the kinematic 

n 

variables [211 ] . We express some scepticism about the possibility of testing LQC effects by using such 
a cosmography. In our opinion, testing the large scale structure (extragalactic matter distribution) 
of the Universe would be more valuable. 



B. Braneworld cosmological models 

The braneworld cosmology is based on the idea that the gravitational field propagates in bulk 
spacetime while our observable universe is only a surface (called brane) embedded in the bigger 
space. This idea leads to the presence of an additional term in the Friedmann equation on the 
brane. (For review of cosmology with extra dimensions see 



22, 



23, 



24 



25]. 



In the present paper we consider a version of higher-dimensional cosmology in the Randall and 
Sundrum framework [gf], [3], which was already tested to some degree with SN la observations 



28, 



In this cosmology, the Friedmann equation assumes the following form 



30, 



3l| 



Tr2 A 4 k 8vr / 4vr \ 2 , C 

where A4 is the 4 dimensional cosmological constant, e = ±1 and C is an integration constant 
whose magnitude as well as sign can depend on the initial conditions. In the present paper we 
take into account the third term of type p 2 . This term arises from the imposition of a junction 
condition for the scale factor of the brane. This condition has simple interpretation - physical 
matter fields are confined to the brane. Both negative and positive e are possible mathematically 



because e corresponds to the metric signature of the extra dimension 

Till now even the sign of e remains an open question. Sahni and Shimmy 32] discussed some 
consequences of a special choice of e = — 1. On should note that possibility that the parameter has 
negative sign is crucial because in this case models with timelike extra dimension can avoid initial 
singularity by bounce [331] . The presence of (negative) p 2 term in the right hand side of ([3]) lads to a 
contracting universe to bounce instead of big-bang curvature singularity: p — > 00, R a b C( iR abcd — > 00. 
If the bounce take place then we have constraint on the value of brane tension a: \a\ > lMeV 
because the bounce takes place at densities greater than during the nucleosynthesis. Note that 
from observational point of view the modification of type (-)p 2 were also recently investigated 



by authors in the context of two-brane model [34]] and this corrections also arises in classically 
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constrained gravity 35|, l36|] . 

The fifth term in ([3]) (called dark radiation) we put here equal to zero. Note that this term 
would decay as rapidly as a -4 and would not modify the early evolution as significantly as the 



third term. However, p 



on nucleosynthesis 



37, 



ease note that non zero dark radiation term can have non-negligible impact 



38] 



C. Non-Riemannian cosmological models 

The cosmological models basing on the so called non-Riemannian extension of general relativ- 
ity theory have been studied for long time (for a review see 4, l^ . I40I]). Trautman introduced 
the first non-Riemannian cosmological models based on Einstein-Cartan theory [4l|. They are a 
modification of general relativity theory by adding the torsion of the spacetime. In such models, a 
consequence of spin and torsion is adding to the Friedmann equation an additional term p 2 oc a~ 6 
The main advantage of such a model is that the problem of initial singularity can be avoided due 
to spin effects. Among the different extensions of General Relativity, the so called metric- affine 
gravity (MAG) is recently being studied. In contrast to Riemann-Cartan theory, the connection is 
not longer metric which implies that covariant derivative of the metric does not vanish. 



Recently some authors 



42, 



2, 



44| have used the magnitude-redshift relation and observations 
of distant supernovae type SNIa to constrain the model parameters. In the case considered the 
modified Friedmann equation assumes the following form 

3 a 2 3 + a 6 ' lj 

where the new constant v can be both negative or positive. The above equation can be rewritten 
to a new form: 

^k,0 + ^m,0 + ^A,o + ^ip,0 = 1> (5) 

where O^o = Vjp, = f^oa" 6 and tp = const is the integration constant. Therefore non- 
Riemannian quantities in this model (torsion and non metricity) will modify very early stages of 
evolution and are negligible at late times. The additional density parameter fl^o is related to the 
non-Riemannian structure. 

Different astronomical observations from Fanaroff-Riley Type lib radio galaxies (FRIIb RG), X 
ray gas mass fraction to SNIa data can be used to constrain the present value of the parameter 
(p.${z = 0) = Slipfi). This gives limits to the possible non-Riemannian structure of the spacetime. 
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The above are main possible interpretations of the presence of a p 2 contribution in the Friedmann 
equation. We must remember that cosmography which bases on the behaviour of null geodesies 
maps the geometry and kinematics of the Universe in terms of H{z) without reference to the 
particular structure of each contribution. Because it measures only average properties of matter 
density it is possible to come up with different forms of contribution leading to the same form 
(— )(1 + z) G [451 ] . It was also showed in 45j that in the generic case, alpha varying models lead to a 
bouncing universe. Interestingly, similar a -6 modifications have also be obtained in the universes 
with varying constants 461 ] . 



All the models under consideration can be represented in terms of density parameters For 
example for the brane models we can use a set of parameters 

o k _-_h_ o P o - _^±_ o _ a p2 

si fc — o tto — -9' ni — o tto ' A — o tto > dr — „tt9) s 'mod 



a 2 H 2 d 2' " m 3 #2' ""A 3jH - 2 ' 3^2' -™ 3H 2 pcr 

where a is the brane tension. Then equation ([3]) can be rewritten to the form [47| 



Note that m and mo d are not independent, i.e. ^modo = ~ o m '° ; ^loops o = —^fh- Moreover 

1 "loons. fl " ' oiir\ 



9 2 

^loops,0 

^ioops,o is fixed from theory if Hq is known. With Hq = 65 km s _1 Mpc -1 and £l m fi — 0.3 it gives 
^ioops,o — 5.24 x 10 122 which gives a contribution to ^ mo d,o — 1-72 x 10 -124 only. It would be 
worth mentioning that there are other sources of correction which act as stiff matter for example 
self-gravitational corrections of an Ads black hole 48]. 



In our paper we consider standard matter with energy density p and pressure p = wp, w = const 
rather than matter in the form of non-minimally coupled to gravity scalar field. Bojowald advocated 
the latter approach where matter emerges from the scalar field [49[]. We assume that standard 
matter satisfies the conservation condition which will determine the dependence of energy density 
on scale factor (or redshift). This means that we are not considering the inverse volume effects 
that would modify the scalar field energy and then conservation conditions. Good news for our 
estimation is that for the case of closed model inverse volume effects are negligible if universe is on 
macroscopic level lo| . 
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FIG. 1: The 68.3% and 95.4% confidence levels (obtained from combined analysis of 
SN+RG+SDSS+CMBR) on the (O m , ,O A , ) plane. 

II. OBSERVATIONAL CONSTRAINTS ON THE FRW BOUNCING MODEL'S PARAM- 
ETERS 



In comparison of the model with observational data we consider two different strategies. First 
of all, because the origin of the p 2 term is not a priori known we fitted a general f2 mo d,o parameter 
from the observations. Accordingly, we obtained that it is too small for detection (probing) by 
astronomical observations (cosmography). In the second approach, we took into account the fact 
that in the LQG, the density parameter f2 mo d,o is fixed. We consequently obtained that the model 
considered is indistinguishable from the ACDM model. It is a consequence of the theory itself 
determining the additional parameter's value to be very small. 

Cosmological models are frequently tested against supernovae observations using the luminosity 
distance d,L of the la supernovae as a function of redshift [1]. With these types of tests for distant 
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SNIa, we can directly observe not the luminosity distance di but their apparent magnitude m and 
redshift z. Taking into account the fact that absolute magnitude M of the supernovae is related 
to its absolute luminosity L, we can obtain the following relation between the distance modulus fi, 
the luminosity distance, the observed magnitude m and the absolute magnitude M 

fi = m-M = 51og 10 d L + 25 = 51og 10 L> L + .M (8) 

where Dl = HqcIl and M = —5 log 10 -^o + 25. The luminosity distance of a supernova is a function 
of redshift and can be computed from the formulae: 

dl( ^ (1+z) i7raK*°^Iw)) (9) 



where 



H\ 2 



— = n m , (i + zf + n fc0 (i + z) 1 + ft modj0 (i + zf + n A)0 , (10) 

Hq) 

^fc,o = — jp and F(x) = (sinh(x), x, sin(x)) for k < 0, k = 0, k > 0, respectively. 

We discuss, for greater generality, non-flat models and keep in mind that for LQG only the 
results obtained with fixed, zero k apply. This happens because LQG models' simple correction 
of the p 2 type appear in spatially flat universes, and assume a more complicated form for k = ±1 



Daly and Djorgovski [5l( (see also [43, |52j, |53|]) suggest including in the analysis not only 



supernovae but also radio galaxies. They pointed out that for tests based on radiogalaxies it is 
useful to apply not the luminosity distance di but the coordinate distance y(z) [3]. The relation 
between the luminosity distance di and the coordinate distance y(z) has the following form 

yy ' c(l + z) y J 

Daly and Djorgovski [55] have compiled a sample comprising the data on y(z) for 157 SNIa in 

the Riess et al. |56l Gold dataset and 20 FRIIb radio galaxies. In our data sets we also include 

□ 

115 SNIa compiled by Astier et al. [571 ]. 

^From the comparison of eq. ()10p and (jllj) it is easy to see that the coordinate distance y(z) 
does not depend on the value of Hq. Unfortunately, we do not observe the coordinate distance y(z) 
of SNIa directly. This distance must be computed from the luminosity distance (or the distance 
modulus n). It is clear that for such a computation a knowledge of the value of Hq is required. For 
both supernovae samples we choose the values of Hq which were used in the original papers. We 
used the distance modulus presented in Ref. js^l. for the calculation of the coordinate distance. 



For each sample we choose the values of Hq appropriate to the data sets. For Riess et al.'s Gold 
sample we have h = 0.646 as the best fitted value and this value is used for calculation of the 
coordinate distance for SNIa belonging to this sample. In turn, the value h = 0.70 was assumed 
in the calculations of the coordinate distance for SNIa belonging to Astier et al.'s sample, because 
the distance moduli [i presented in Ref. 57|, Tab. 8] were calculated with such an arbitrary value 



of h = 0.70. 

The error of the coordinate distance can be computed as 

\ 2 / 



10" 



c(l + z) 10 5 



Hq In 10 



(12) 



where Cj(/ij) denotes the statistical error of the distance modulus determination (note that for 
Astier et al.'s sample the intrinsic dispersion was also included) and ct(Hq) = 0.8 km/s Mpc 
denotes the error in Hq measurements. 

With the use of coordinate distance y(z) as our basic quantity it is easy to include in our analyse 
two additional constraints which do not depend on the value of Hq either. These constraints are 
obtained from extragalactic analysis. First, we have the baryon oscillation peaks (BOP) detected 



in the Sloan Digital Sky Survey (SDSS) Luminosity Red Galaxies 
A 



A 



\/^m,0 




They found that value of 



(13) 



Sfc,0| \ ' Jo 

(where E(z) = H(z)/Hq and z\ = 0.35) is equal to A = 0.469 ± 0.017. The quoted uncertainty 
corresponds to one standard deviation, where a Gaussian probability distribution has been assumed. 

The second constraint which we include in our analysis is the so called (CMBR) "shift param- 
eter" 




dz 



where Rq = 1.716 ± 0.062 



59] 



In our combined analysis, we can obtain a best fit model by minimizing the pseudo-x 2 merit 
function [6C 



2 _ 2 | 2 | 2 

X — XSN+RG ~r XSDSS + XCMBR 



E 



„ obs 



y th\2 /^mod _ o.469 \ 2 /R mod - 1.716 
+ V 0.017 J + V 0.062 



(14) 
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where ^4 mod and R mod denote the values of A and R obtained for a particular set of the model 
parameter. For Astier et al.'s SNIa sample [57I] an additional error in z measurements was taken 
into account. Here di(yi) denotes the statistical error (including the error in z measurements) of 
the coordinate distance determination. 

We can obtain constraints for the cosmological parameters by minimizing the following likelihood 
function C oc exp(— \ 2 /2). One should note that when we are interested in constraining a particular 
model parameter, the likelihood function marginalized over the remaining parameters of the model 
should be considered [60(]. 

Our results are presented in Table HI Table [H] and Fig. [TJ Table U refers to the minimum x 2 
method, whereas Table HI1 shows the results from the marginalized likelihood analysis. In Tabic UTTl 
Table HV] we repeated our analysis with the prior f2 mo d,o — 0- 

^From our combined analysis (SN+RG+SDSS+CMBR), we obtain as the best fit a flat (or 
nearly flat universe) with O mj o — 0.3, and Oa.o — 0.7. For the (l + z) 6 term we obtain the stringent 
bound fi mo d,o ^ (—0.26 x 10~ 9 ,0.31 x 10~ 4 ) at the 95% confidence level. These results mean 
that the positive value of f2 m od,o is preferred (f2 m od,o > 0), however small negative contribution 
of (1 + z) e type is also available. Our results shows that in the present epoch contribution of the 
dark radiation, if it exists, is small and gives only small corrections to the ACDM model in the low 
redshift. 

The above analysis was performed for any spatial curvature. However, please note that it 
holds for non-flat models only if the Friedmann equation does not change in any other way after 
introducing the p 2 term. Unfortunately, LQG models with non-zero curvature were obtained 
for different setup of matter components 



13, 



50|]. This means that for LQG only the flat case 



consideration apply. 

The general basic formula for H(z) for universe with Qk,o > rewritten in terms of dimensionless 
density parameters f2, = is 



(i 



~ET = ( fi m,o(l + 20 +fifc,o(l + z) +^A,o) 1 p; 3 , (15) 

V-TCO/ V "loops.O "loops^ / 



a 2 

(where S] mo do = — n m '° ). One should note that in the case considered A also enters in non- 

^loops,0 

standard fashion. It means that in this case we do not have the simple reduction to one model 
with different interpretations of S^ moc j, and the same Friedmann equation. Then, we are no longer 
comparing simple (1 + z) 6 modifications, but rather different curved models. This case (under the 
prior ftkfi > 0) are analysed and presented separately. 
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When one wants to compare the results obtained for different models, an interesting question 
is: how significant is the improvement of the fit due to the new model? To answer this question 
one can use the information criteria. The most popular information criteria used in everyday 
statistical practise are the Akaike information criterion (AIC) 



611 ] and the Bayesian information 



criterion (BIC) 62|. These criteria can be used also for selection of model parameters by providing 
the preferred data fit. 

Usually, incorporating new parameters increases the quality of the fit. The question is, if it 
increases it significantly enough. The information criteria put a threshold which must be exceeded 
in order to assert an additional parameter to be important in explanation of a phenomenon. The 
discussion of how high this threshold should be caused the appearing of many different criteria. For 
us, it suffices to check whether the AIC and BIC provide sufficient arguments for incorporation of 
the new parameters. The power of u sing information criteria of model selection was demonstrated 



by Liddle 



631 ] and Parkinson et al. 



6j|]. Please note that in any case some future observational 



data may give arguments in favour of additional parameters. Also, theoretical considerations could 
lead to inclusion of new parameters, even it is not necessary from the point of view of the present 
observations. 

The AIC is defined in the following way 61] 



AIC 



-21n£ + 2d, 



(16) 



where C is the maximum likelihood and d is a number of free model parameters. The best model 
with a parameter set providing the preferred fit to the data is the one that minimizes the AIC. 
It is interesting that the AIC also arises from an approximate minimization of the Kulbak-Leibler 
information entropy 



65]. 



The BIC introduced by Schwarz 



621 ] is defined as 



BIC 



-2\nC + d\nN, 



(17) 



where N is the number of data points used in the fit. Comparing these criteria, one should note that 
the AIC tends to favour models with a large number of parameters unlike the BIC, which penalizes 
new parameters more strongly. It is the reason why the BIC provides a more useful approximation 



to the full statistical analysis in the case of no priors on the set of model parameters |64l ]. This 
means that while the AIC is useful in obtaining upper limit to the number of parameters which 
should be incorporated to the model, the BIC is more conclusive. One should note that only the 
relative value between the BIC of different models has statistical significance. The difference of 2 is 
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treated as a positive evidence (and 6 as a strong evidence) against the model with the larger value 



of the BIC 



66 



671] - If we do not find any positive evidence from information criteria the models 
are treated as a identical and eventually additional parameters are treated as not significant. The 
using of the BIC seems to be especially suitable whenever the complexity of reference does not 
increase with the size of data set. The problem of classification of the cosmological models on the 
light of information criteria on the base of the astronomical data was discussed in our previous 



papers 



68, 



69 



70, 



71 



Our results are presented in Table PVT. Please note that the results of statistical analysis for LQG 
models with non-zero curvature (under prior Qj^o > 0) are presented as a separate case. It is clear 
that in the light of information criteria, the (1 + z) 6 term does not increase the fit significantly. It 
confirms, what could be expected, that this term, if it exist, is small in the present epoch 

In the Bayesian framework, the quality of models can be compared with help of evidence 



We can define the a posteriori odds for two models - Mi and Mj - the so called Bayes factor B 



66|,l24]. 
3- 



If a priori we do not favor any model it reduces to the evidence ratio. Schwarz |62| showed that for 
observations coming from a linear exponential distribution family in the asymptotic approximation 
N — > oo the logarithm of evidence is given by 

\nE = lnC-^\nN + 0{l). (18) 

It is easy to show that in this case we have a simple relation between the Bayes factor and the BIC 

2 In Bij = -{Bid ~ BICj) (19) 

If Bij is greater than 3 it is considered a positive evidence in favor of Mj model, while B^ > 20 



72]. We present our 



gives strong, and B^ > 150 very strong evidence in favor of model Mj 
results in Table IVI1 In all cases we obtain positive evidence in favour of the ACDM model over 
the bouncing cosmology with the term p 2 . Only in the case when additional parameter of the 
theory (responsible for the term (1 + z) & ) is fixed the obtained model is indistinguishable from the 
ACDM model. It supports the results obtained with help of the AIC and BIC. This result is valid 
both with and without the prior of r2 mo d,o — 0- One should note that because of a non-Gaussian 
distribution of the a posteriori PDF function for £l mo d,o the results obtain with help of Bayes factor 
(especially with the priors fi mo d,o — 0) should be treated with caution as only additional support 
for results obtain with the AIC and BIC. 

Please also note that if ^ m od,o < 0, then we obtain a bouncing scenario [3, H, H] instead of 
a big bang. For O m) o = 0.3, f2 m od,o = —0.26 x 10~ 9 and h = 0.65 bounces (H 2 = 0) appear for 
z ~ 260. In this case, the BBN epoch never occurs and all BBN predictions would be lost. 
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III. LOCATION OF CMB PEAKS AND BBN IN THE MAG MODEL 



The results obtain in the previous section lead to the conclusion that we should obtain stronger 
constrains for model parameter especially for fJ mo d.o- One should note that for doing so, it is useful 
to analyse the location of the first peak in the CMB power spectrum and the predictions of the 
BBN. Stronger constraints for model parameters in the MAG model was obtained in Ref. [441 ] . 

The idea of testing models using the location of the first peak in the CMB power spectrum 
is based on the fact that the hotter and colder spots in the CMB can be interpreted as acoustic 
oscillation in the primeval plasma during the last scattering. Peaks in the power spectrum corre- 
spond to maximum density of the wave. In the Legendre multipole space these peaks correspond 
to the angle subtended by the sound horizon at the last scattering. Further peaks correspond to 
higher harmonics of the principal oscillations. The locations of these peaks are very sensitive to 
the variations in the model parameters. Therefore, the position of the first peak can be used as 
another way to constrain cosmological models. 

In the MAG model, assuming f2 mi o = 0.3 and h = 0.72 for the standard ACDM universe, the 
correct positions of the first peak was obtained in [44j] as t\ = 220. ^From the SNIa data analysis, 
it was found that the Hubble constant has lower value. Assuming that Hq = 65 km/s Mpc (or 
h = 0.65) and consider the standard ACDM model, with r2 mj o = 0.3, one gets l\ = 225 44]. 

Some discrepancy between the observational and theoretical results was found in this case. Now 
it was interesting to check whether the presence of the fictitious fluid fi^o changes the locations of 
the peaks. If we choose the Hq = 65 km/s Mpc then agreement with the observation of the location 
of the first peak could be obtained for three non-zero values of the parameter O^o- Two positive 
and one negative values of this parameter for which the MAG model is admissible are 3 x 10 , 
7 x 10~ 14 and -1.4 x 10~ 10 . 

One should note that agreement with prediction of big-bang nucleosynthesis (BBN) is also 
crucial for testing of the model. Of course the big-bang nucleosynthesis is a very well tested area of 
cosmology and does not allow for any significant deviation from the standard expansion law apart 
from very early times (i.e., before the onset of BBN). The predictions of standard BBN are in good 
agreement with observations of the abundance of light elements. Therefore, all nonstandard terms 
added to the Friedmann equation should give only negligible small modifications during the BBN 
epoch to leave the nucleosynthesis process unchanged. 

In our opinion the consistency with BBN is a crucial issue in the models where the nonstandard 
term o -6 is added in the Friedmann equation. It is clear that such a term has either accelerated 
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(^mod,o > 0) or decelerated (r2 mo d,o < 0) impact on the Universe expansion. Going backwards 
in time this term would become dominant at some redshift. If it had happened before the BBN 
epoch, the radiation domination would have never occurred and all the BBN predictions would be 
lost. 

If we assume that the BBN result are preserved in our model, we obtain another constraint on 
the amount of Ow, o- Let us assume that the model modification is negligibly small during the BBN 
epoch and the nucleosynthesis process is unchanged. It means that the contribution of the term 
^■mod cannot dominate over the radiation term fi r; o ~ 10~ 4 before the beginning of BBN [z — 10 8 ) 

fWi,o(i + ^) 6 < n rfi (i + zf |^, i0 | < io- 20 . 

The values of fi m odo °c 10~ 3 obtained as best fits in the SNIa data analysis as well as the 
smallest nonzero value of fi^o = 7 x 10~ 14 calculated in the CMB analysis are unrealistic in the 
light of the above result. If we take into consideration the maximum likelihood analysis of SNIa 
data we have the possibility that the value of Q mo d,o is lower than 1 10 20 1 in the 2a confidence 
interval. 

One could argue that bounds from cosmography and CMB are weaker than those from nucle- 
osynthesis. However, they are model independent. Of course, nucleosynthesis is a well tested area 
in cosmology [79| but it is described rather in the terms of standard physics without loops correc- 
tions. The first s tep toward the description of Big-Bang nucleosynthesis are given in the paper by 
Bojowald et al. [801 ] . The authors demonstrated that several correction to the equation of state 
parameter can arise from classical and quantum physics, for example loops quantum gravity allows 
one to compute quantum gravity corrections for Maxwell and Dirac field. So we cannot a priori 
assumed that the corresponding gravity corrections are negligible during the Big-Bang nucleosyn- 
thesis. In the present authors' opinion, the significance of bounds obtained from cosmography and 
CMB is that they are independent from the nucleosynthesis bounds and are related to different 
cosmological epochs. 

IV. CONCLUSION 

In the paper we have studied observational constraints on the FRW models with p 2 modifica- 
tions, and obtained stronger limits on the magnitude of the term fl mo <i (scaling like (1 + z) G ). We 
pointed out that astronomical observations allow us to test the total contributions of a fluid scaling 
like (1 + z) 6 but we cannot separate particular terms of a negative sign. It is a simple consequence 



14 



TABLE I: Results of the statistical analysis of our model with (1 + z) 6 like term obtained from \ 2 best fit. 
The upper section of the table represents the constraint Qk,o = (flat model). 

sample Clk,o ^m,o ^mod,o ^a,o X 2 

SN - 0.23 0.020 0.757 295.7 

SN+RG - 0.26 0.013 0.727 319.5 

SN+RG+SDSS - 0.28 0.009 0.711 319.6 

SN+RG+SDSS+CMBR - 0.30 0.4 x 10~ 8 0.700 322.4 

SN 0.34 0.00 0.044 0.646 295.6 

SN+RG 0.29 0.05 0.034 0.626 319.5 

SN+RG+SDSS -0.03 0.28 0.011 0.739 319.5 

SN+RG+SDSS+CMBR 0.04 0.30 0.279 x 10~ 6 0.660 321.4 



TABLE II: Results of the statistical analysis of our model with (1 + z) e term. The values of the model 
parameters are obtained from marginalized likelihood analysis. We present maximum likelihood value with 
68.3% confidence ranges. The upper section of the table represents the constraint Clkfi = (flat model). 

sample fifc,o fi m ,o ^mod,o ^A,o 

SN - 0.22±g : g5 0.0211°;°?° 0.75±°;°! 

SN+RG - 0.27l°;°| 0.0131°;°^ 0.731°,$ 

SN+RG+SDSS - 0.281°;°^ 0.0091°,;°,°^ 0.7lt°,;jg 

SN+RG+SDSS+CMBR - 0.30l°;°? (0.279 x lO^ 8 ) :^ 1 ^^ 0.701°;°^ 

SN 0.10l°j° 0.00i°;g* 0.024t°;°£ 0.77±°;i| 



SN+RG 0.14l°;J 8 O.OOlg;" 0.0191°;°^ 0.73t°;i$ 

SN+RG+SDSS -0.03±g;?g 0.28t° o f 2 0.01ll°;°^ 0.74±g;jg 

SN+RG+SDSS+CMBR 0.0ll°;°4 0.30l°;°? (0.36 x lO" 6 )!^*^ 6 0.69t°;°^ 



of the fact that cosmography measures only the "average" density of matter. However, we showed 
that some stringent bounds on the value of this total contribution can be given. 

There are several interpretations of the presence of the (1 + z) 6 term in the Friedmann equation: 
1) the brane theory; 2) the non-Riemannian theory of gravity; and 3) the loop quantum cosmology. 
All models give rise to a (1 + z) 6 correction which effects are important in the very early universe 
and become unimportant in later evolution. The loop quantum cosmology is more fundamental 
theory than the classical theory. We are looking for the quantum cosmology predictions which are 
a priori in agreement with the classical picture of the present Universe (we showed the quantum 
loop effect are negligible at the present epoch). 

Note that if we consider the phantom p 2 modification, it is also important in the future evolution 
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of the universe (Big-Rip singularities) 



We used Bayesian methods of model selection to answer the question: which cosmological model 
- with initial singularity or with bounce - is promoted by observational data? We have shown that 
models with a singularity are a "more economical" choice for the Universe but bouncing cosmology 
cannot be ruled out even on the la level. However, because the the posteriori probability function 
for f2 mo d,o is strongly non-Gaussian, probability of bounce, i.e. -P(f2 m od,o < 0), is less than 1%. 
Our observational analysis clearly reflects the important role of independent observational data 
which enables us to refine the analysis of model parameters. 

The analysis of SNIa data as well as both SNIa and FRIIb radio galaxies (with and without 
priors coming from baryon oscillation peaks and CMBR "shift parameter") shows that the values 
of x 2 statistics are lower for model with (1 + z) 6 like term, than for the ACDM model. On the 
other hand, information criteria show that including such a term does not increase the quality of 
the fit significantly; or, alternatively, that the quality of the available data is not good enough for 
fitting this new term. BIC even favours the ACDM model over our model of bounce, although this 
preference is weak. These results lie in agreement with the fact that (1 + z) 6 term is not significant 
in the present epoch of the Universe. Moreover, if the additional parameter of the theory f2 mo d o is 
fixed, like for the LQG model, than we obtained model is indistinguishable from the ACDM model 
in the Bayesian framework. 

The combined analysis of SNIa data and FRIIb radio galaxies using baryon oscillation peaks 
and CMBR "shift parameter" gives rise to a concordance universe model which is almost flat with 
O m o — 0.3. ^From the above mentioned combined analysis, we obtain the following constraint for 
the term which scales like (1 + zf: ft mod ,o e (-0.26 x 10~ 9 ;0.31 x 10~ 4 ) 

We confirm Bojowald's assertion that effects of quantum gravity can by potentially tested, but 
the required bounds lie beyond the possibilities of high precision cosmology. What we obtain from 
cosmography is very weak and other tests may by useful. We have shown that the analysed scenario 
is compatible with the most recent low redshift observations of SN la, which are independent of 
physical processes in the early universe. We find the other limits on the value of fl^o from 
measurements of CMB anisotropics and BBN. We obtain the strongest limits in this case, namely 
^i/>,o — 10 -20 from BBN. Of course BBN as well as CMB are a very well tested areas of cosmology 
which do not allow for significant and substantial changes. Still, one must remember that, although 
consistency with BBN and CMB is a crucial issue, in such an approach we a priori assume that 
brane models (for example) do not change the physics of the pre-recombination epochs. 

The analysis was performed for flat and non-flat models with the tacit assumption that the 
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TABLE III: Results of the statistical analysis of the bouncing model with the negative (1 + z) 6 term obtained 
from x 2 best fit. The upper section of the table represents the constraint O^o = (flat model). 



Sample fij^o ^ m ,o ^mod,o ^a,o X 

SN - 0.31 0.000 0.69 297.5 

SN+RG - 0.32 0.000 0.68 320.4 

SN+RG+SDSS - 0.30 0.000 0.70 322.4 

SN+RG+SDSS+CMBR - 0.30 0.000 0.70 322.5 

SN -0.28 0.43 0.000 0.85 296.0 

SN+RG -0.18 0.39 0.000 0.80 319.5 

SN+RG+SDSS 0.06 0.29 0.000 0.65 321.1 

SN+RG+SDSS+CMBR 0.00 0.30 0.000 0.70 322.5 

curvature does not change the physics in any other way than introducing an appropriate term in 
the Friedmann equation. Since this is not the case for LQG, only the flat case considerations giv e 
bounds for the parameter Sli oops . Although curved LQG models have been introduced [l3|, l50l |. 
they have significantly different matter components, making it impossible to compare them in our 
scheme. 

Our general conclusion is that while cosmography can be generally used to test p 2 type contri- 
bution, such correction turns out to be very small. Moreover, because we use H(z) function which 
probes only the average density, it is not possible to separate effects of LQG from other effects 
scaling like (1 + z) 6 , for example effects of the brane. On should note that in LQG theory f2i O ops,0 is 
fixed and gives an odd contribution to f2 mo d,o — 1-7 x 10~ 124 only, which is far below the possibility 
to test by present cosmography. It mean that any positive evidence for non zero (1 + z) 6 term 
eventually obtained from cosmography can not be connected with LQG. 
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TABLE IV: Results of the statistical analysis of the bouncing model with the negative (1 + z) e term. The 
values of the model parameters are obtained from the marginalized likelihood analysis. We present the 
maximum likelihood values with 68.3% confidence ranges. The upper section of the table represents the 
constraint Sl^.o = (flat model). 



Sample 



fe,0 



a, 



^A,0 



SN 




- 


n oq+0.02 
u -°°-0.02 


0.000 


-0.006 


0.67 


SN+RG 






r, qq + 0.03 

U.OO_ 02 


0.000 


-0.006 


0.67 


SN+RGH 


-SDSS 




n 31 + 002 

u - ol -0.02 


0.000 


-0.002 


0.69 


SN+RG+SDSS+CMBR 




30+ 02 


0.00000_ 


0.22X10- 9 


0.70 


SN 




-0 42+ - 24 

u -* z -0.22 


n rr,+0.09 

u - oz -0.09 


0.000 


-0.011 


0.91 


SN+RG 




-0.36+o.jo 


n 49+Q-io 
u -^ y -o.n 


0.000 


-0.011 


0.85 


SN+RGH 


-SDSS 


o.o8±85g 


u - zy -0.02 


0.000 


-0.004 


0.63 


SN+RGH 


-SDSS+CMBR 


o.oolg : g 


n sn+o- 02 


0.00000_ 


0.21X10- 9 


0.70 



+0.02 
-0.02 

+0.02 
-0.03 

+0.02 
-0.02 

+0.01 



+0.11 
-0.15 



+0.12 
-0.13 



+0.05 
-0.06 



+0.02 
-0.02 



TABLE V: The values of AIC and BIC for the ACDM model and Bouncing Cosmology model (with the 
term (1 + z) e ) without and with priors r2 mo d,o < 0. Separately the LQG model with fixed f2 mo d,o (^fc,o > 0) 
is considered. The upper section of the table represents the constraint fij^o = (flat model). 

ACDM ABCDM ABCDM^o < 0) LQG 



sample 


AIC BIC 


AIC BIC 


AIC BIC 


AIC BIC 


SN 
SN+RG 
SN+RG+SDSS 
SN+RG+SDSS+CMBR 


299.5 303.1 
322.4 326.1 

324.4 328.1 

324.5 328.2 


299.7 306.9 

323.5 330.9 

323.6 332.0 
326.4 333.8 


301.5 308.7 
324.4 331.8 

326.4 333.8 

326.5 333.9 


299.5 303.1 
322.4 326.1 

324.4 328.1 

324.5 328.2 


SN 
SN+RG 
SN+RG+SDSS 
SN+RG+SDSS+CMBR 


300.0 307.2 
323.5 330.9 

325.1 332.5 
326.5 333.9 


301.6 312.4 
325.5 336.5 
325.5 336.5 
327.4 338.4 


302.0 312.8 
325.5 336.5 

327.1 338.1 
328.5 339.5 


301.5 308.7 

324.4 331.8 
325.1 332.5 

326.5 333.9 



Appendix. Basics of Loop Quantum Cosmology. 



In this appendix we give some selected information about Loop Quantum Gravity (LQG) con- 
nected to the subject of this work. Our main challenge is to show how to obtain equation (1) 



82(|). We concentrate rather on the main steps 



considered in this paper (see also a nice summary 
and for detailed calculations we direct to the references. 

Loop Quantum Gravity describe the gravitational field as a SU(2) non-Abelian gauge field using 
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TABLE VI: The values of Bayes factor for models: 1) ACDM model, 2) Bouncing Cosmology model (with 
the term (1 + z) 6 ) and 3) Bouncing Cosmology model (with (1 + z) 6 term) with priors f2 moc j,o < 0, 4) LQG 
model with fixed r2 mo d,o (^fc.o > 0)- The upper section of the table represents the constraint Clk.o — (flat 
model). 



sample B 12 B 13 B 14 B 23 B i2 B 43 

SN 6.69 16.44 1.00 2.46 6.69 16.44 

SN+RG 11.02 17.29 1.00 1.57 11.02 17.29 

SN+RG+SDSS 7.03 17.29 1.00 2.46 7.03 17.29 
SN+RG+SDSS+CMBR 16.44 17.29 1.00 1.05 16.44 17.29 
SN 13.46 16.44 2.12 1.22 6.36 7.77 

SN+RG 16.44 17.29 1.57 1.00 10.49 11.02 

SN+RG+SDSS 7.39 17.29 1.00 2.23 7.39 17.29 
SN+RG+SDSS+CMBR 9.49 17.29 1.00 1.73 9.49 17.29 



background independent methods. The canonical fields are so called Ashtekar variables (^4, E) [83 ] 
which take value in su(2) and su(2)* algebras respectively. These variables are analogues of the 
four potential and electric field in electrodynamics. The Ashtekar variables are strictly connected 
with triad representation. However in LQG gauge fields describe only spatial part X when time 
is treated separately. To quantise this theory in the background independent way ones introduce 
holonomies of connection A 

h a [A]=Vexp / A where 1-form A = r i A i a dx a (20) 

J a 

and conjugated fluxes 

F i s [E]= I dF l where 2-form dFi = e abc E?dx b A dx c . (21) 
Js 

which are background independent observables. In equation (|20p we have introduced t% = — 
where <7j are Pauli matrices. The quantisation of theory lead to important result that volumes and 
areas have discrete spectrum. For further calculations it will be important to notice that there 

n n 

exists a minimal area A = 2^/3^7^ [8J|. For introductory review on loop quantisation see [851 ]. 
Dynamics of theory is contained in the scalar constraint 

H ° = lib L 63xNix) ^m \- £ " tFi " 2(1 + 72)Kf * A <i] (22) 

where F is a field strength F = dA + \ [A, A] and K is an extrinsic curvature. Constant 7 in 
equation (|22|) is so called Barbero-Immirzi parameter, 7 = ln2/(7T\/3). 
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We want to show now how to apply LQG to case of the flat FRW model considered in this paper. 
The results presented below base on papers where reader can find detailed calculations 

and analysis. The FRW k = spacetime metric can be written as 

ds 2 = -N 2 {x)dt 2 + q ab dx a dx b (23) 

where N(x) is the lapse function (here we choose gauge N(x) = 1) and the spatial part of the 
metric is expressed as 

q ah = Syu&i = ^{tfq^u^ul (24) 



In this expression °q ab is fiducial metric and °oj l a are co-triads dual to the triads °ef, uj l (°ej) = 8\ 
where °lo % = °LO % a dx a and °6i = °efd a . In the case considered the Ashtekar variables are 

A = Y + 1 K = cV^ 1/S °J a Tidx a , (25) 
E = V|detg|e = pV~ 2/3 ^q °e?r<a« (26) 

where Vq is volume of fiducial cell and Y is spin connection. The parameters (c, p) can be since now 
considered as canonical variables with fundamental Poisson bracket {c, p} = 8irGj/3. In quantum 
theory p is replaced by operator p which acting on the eigenvector \fj) gives 

AmHm^Im). (27) 

where \x £ M and eigenvectors fulfils relation of orthogonality (/Jti\fij) = S^.^. There is no well 
defined operator of variable c, instead of this there exists another fundamental operator defined as 

exp^l/x) = \fjL + X). (28) 
With use of definition (|20j) we can calculate holonomy for connection (I25p in particular direction 

^=Icos(^)+r,sin(^). (29) 

Holonomy is well defined operator which acting on vector gives 

^V) = \ (lM + A) + l/x - A)) + ~ (\n + A) - |/i - A)) Tl (30) 

where we have used definition (I28p . lYrom particular holonomies (I29D we can construct holonomy 
along the closed curve a = EL-. This holonomy can be written as 

h M = ftO.) h M^)-^M-l = n + fV^F k ab r k °ere) + 0(„ 3 ) (31) 
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and inverting this equation we obtain expression for field strength 



tr 



F 



Tk [h 



ab 



o to J 



M 2 ^o 2/3 P 2 < 



sin 2 (/ic) j 



(32) 



In the limit /i -t we retrieve classical formula. However this limit does not exist in quantum 
theory because of area gap A. LQG tell us that we should stop shrinking the loop /i — > p, when 
physical area pfl 2 = A. Then inserting expressions for field strength (|32l) and for A and K to (|22|) 
we obtain the effective Hamiltonian with holonomy corrections 



8irGr 



■y/P 



sin (/tc) 



(33) 



where we have also added matter part. Additional feature of theory is Hamiltonian constraint 
H e Q = which implies 



1 

7 2 p 



sin (/ic) 



8vrG 



Now we can use the Hamilton equation to determinate dynamics of the canonical variable p 

2 fp 

p = {p, H cS } = sin (pc) cos (/2c) . 

7 M 



(34) 



(35) 



Combining equations (|34|) and (|35|) we finally obtain the modified Friedmann equation with holon- 
omy correction 



8ttG 



H 2 

where we have introduced critical energy density 

3 3 

Pc 



Pm 1 



V3 



8vrG7 2 /2 2 p 8vrG 7 2 A 16^1^' 
For further details of Loop Quantum Cosmology we recommend the introductory review 
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